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Abstract 

In this paper we study the effect of stochastic perturbations on a common type of moving 
boundary value PDE's which endorse Stefan boundary conditions, or Stefan problems, and 
show the existence and uniqueness of the solutions to a number of stochastic equations of 
this kind. Moreover we also derive the space and time regularities of the solutions and the 
associated boundaries via Kolmogorov's Continuity Theorem in an appropriately defined 
normcd space. 

The paper first conveys our previous results where randomness is smoothly correlated in 
space and Brownian in time, then introduces a new methodology that enables us to prove the 
existence and uniqueness of the solution to the standard heat equation driven by a space-time 
Brownian sheet as well as its boundary regularity, and finally extends it to the stochastic 
moving boundary partial PDE's driven by the same type of randomness. 

1 Introduction 

An important type of problems in the theory of partial differential equations (PDE's) is the 
moving boundary value problems. In this paper we study the effect of stochastic perturbations 
on such type of problems with Stefan boundary conditions, or Stefan problems, which have 
various applications in physics, engineering, and finance, and show the existence and uniqueness 
of the solutions to a number of stochastic equations of such kind. We also obtain the space and 
time regularities of the solutions and their associated boundaries via Kolmogorov's Continuity 
Theorem in a defined normed space. 

The paper first conveys our previous results where randomness is smoothly correlated in 
space and Brownian in time, which has the convenience that the boundary regularity is natu- 
rally satisfied, then introduces a new methodology that enables us to prove the existence and 
uniqueness of the solution to the standard heat equation driven by a space-time Brownian sheet 
as well as its boundary regularity by defining an appropriate normed space whose norm correctly 
curvatures both the decay of the solution and its regularity at the boundary, and finally extends 
it to the stochastic moving boundary partial PDE's driven by the same type of randomness, 
where the new normed space is defined essentially in the same manner as enlightened by the 
previous result. 

1.1 Background 

A moving boundary PDE of u{t, x) describes the behavior of a system that consists of two 
phases, as illustrated in Figure [l} where f3{t) is a moving boundary which is part of the solution 



m 

Figure 1: An Illustration of Moving Boundary PDE's 

and must be solved simultaneously with u{t, x). As can be seen in Figure [l| in the region to the 
left of the moving boundary (namely, the set {{t, x) £ [0, oo) x M : x < /5(t)}) u is constantly set 
to 0; on the right side of the boundary (the set {{t, x) £ [0, oo) x M : x > f3{t)}) u is described by 
a PDE of the general form Lu + 6 = where L is a predefined second-order differential operator. 

For a moving boundary PDE problem, in addition to the regular boundary condition such as 
the Dirichlet condition u{t, (3{t)) = 0, there is always an extra boundary condition that describes 
the dynamics at the moving boundary, for instance, the Stefan boundary condition 

rill 

-{t,m+) = pm- (1) 

The type of moving boundary PDE's we consider throughout the thesis is the Stefan prob- 
lems, where L is a heat or parabolic operator (for instance, L := —d/dt + d'^/dx'^) with the 
Stefan boundary condition. Such type of problems has a variety of applications. For instance, 
in physics, they model the phenomena such as ice melting with the Stefan condition describing 
the heat balance at the interface (the moving boundary, see [2]); in finance they model the 
valuation of American options with the PDE derived from the Black-Scholes formula and the 
moving boundary describing the early exercise price boundary (see Lemma 7.8, Chapter 2, [4]). 

1.2 Motivation and Results 

In the mathematics of this thesis we are interested in the stochastic versions of the Stefan prob- 
lems, namely, b{t, x) is a formal notation about the stochastic addition (noise). In general b{t, x) 
is a 2-dimensional distribution and therefore we work within the framework of the stochastic 
PDE theory by Walsh in [7], and based on the weak formulations of the equations and their 
equivalent evolution equations. 

When b{t, x) is multiplicative, namely, b{t, x) = u{t, x)W{t, x) where W is the noise (formal), 
[5] proves the existence and uniqueness of the solution when W = W{t) is a distribution (Brown- 
ian) only in time and is constant in space. Also, in |6] we proved the existence and uniqueness of 
the solution when W{t, x) is a distribution (Brownian) only in time and is smoothly correlated 
("colored") in space, which is quickly reviewed in Section [2] 
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However, when W{t, x) is a distribution (Brownian) in both space and time, we face a number 
of novel challenges that are beyond the scope of current literature: 

(1) The spatial derivatives of the solution may not exist (see |7]), which means the techniques 
we used in [6j based on ff-norms may not be available; more importantly, the Stefan 
condition ([T]) involves the spatial derivative of the solution at the boundary, and therefore 
we shall show the existence of such a derivative simultaneously with the existence and 
uniqueness of the solution. 

(2) The stochastic perturbation given by b{t, x) is no longer spatially Lipschitz as in [5] and [6], 
which means in order to control the boundary shift effect in the Ito integrals and the 
nonlinear drift term we may need additional spatial regularities on b{t, x) in addition to 
just being multiplicative as in [5] and [6]; in other words, although the perturbation b 
vanishes when u = 0, u itself may not provide sufficient spatial smoothness to control the 
Ito integrals when the boundary shifts. 

To tackle (1) alone, in Section [s] we first study the stochastic heat equation driven by a 
multiplicative space-time Brownian noise, namely, with W a standard 2-dimensional Brownian 
sheet defined in [7], 

du d'^u d^W 

u(0, x) = uo{x), 

u{t,0) = 0,Vx < 0,t G [o,r]. 

From such study we obtain the conditions and form of norms under which a boundary derivative 
exists, and more importantly, develop the essential techniques to calculate the sample- wise space 
and time modulus of continuity of the solution by means of Kolmogorov's Continuity Theorem 
(see [3j), which is critical to evaluate the regularity needed to control the boundary shift effect. 

Next, in Section [4| first with the same multiplicative space-time Brownian noise (that is, 
b(t, x) = uWtx,) we preliminarily calculate the effect of a boundary shift on the Ito integral using 
the techniques and results obtained in the previous section, and find that it may be difficult 
to obtain the wanted control on iteration. Therefore we make the change in the stochastic 
perturbation b{t, x) such that it is even smoother in space and in turn study the following 
stochastic Stefan problem driven by a scaled space-time Brownian noise: 

du d'^u d^W 

u{t,x) = 0,yx< l3{t), 
u{0,x) = uo{x), 

A/ N -, U(t, X) ^ 

p${t)= lim -^^,P-a.s. 

x\(S{t) X 

where b{t,x) = a{x — li{t))Wtx and (t(-) is a function that satisfies certain regularity conditions 
which are sufficient to tackle (2). We show the existence and uniqueness of the solution to 
the above equation and additionally the regularity of the boundary using results from this and 
previous sections. This result also serves as the mathematical foundation for the modeling of 
the dynamics of limit orders. 
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2 A Stochastic Stefan Problem with Spatially Colored Noise 

In this section we quickly review our work in [6j by giving the main theorems and lemmas 
without proofs. We studied a stochastic Stefan problem of u(t, x) driven by a multipHcative 
noise u{t, x)d^tix), where is a noise that is Brownian in time and smoothly correlated (or 
"colored") in space. Specifically, fix a probability space {Q, 3", P), and suppose : Q x M x M — )■ 
M is the standard 2-dimensional Brownian sheet. Suppose also that r/ : M — >■ M is C°° and 
= 1. Then for t > 0, x G M, define 

pt poo 

Ctix) ■■= / r]{x-y)W{dyds). 
Jo Jo 

Then we have the following problem and results. 
2.1 Problem and Main Theorem 

We showed the existence and uniqueness of the solution u{t, x) to the following formal equation: 

~qI = + u{t,x)dit{x),yx > Pit), 

u{t,x) = 0yx<(3{t), 
u{0,x) = uo{x), 

pp(t)= lim ^^i^,P-a.s. 

x\l3{t) X 

for < t < r where r is some well defined stopping time. 

Since 1^ is a distribution, equation ^ is in fact formal, and we need to work on the weak 
definition (Definition 3.2 in |6]) and its equivalent evolution equation (Equation (15) or Lemma 
3.4 in [6j). The main theorem is 

Theorem 2.1 The solution u{t,x),(3{t) to Q exists and is unique for for < t < t := 
liiQL^^T^ where := inf{t G M+ U {0} : |/3(i)| > L}, and u{t,x) := u{t,x + (3{t)) satis- 
fies 



oo 



u{t,x)= p{t,x,y)uoiy)dy 





ft fOO 

+ ${s)q{t - s,x,y)u{s,y)dyds (3) 

Jo Jo 

ft poo 

+ 11 pit - s,x,y)u{s,y)d^tix)ds, 
Jo Jo 



where p, q are standard kernels defined in f^. 



2.2 A Transformation 

We make the natural transformation uit,x) := uit,x + /3(t)), which transforms the original 
weak definition and its equivalent evolution equation to a nonlinear PDE in the fixed domain 
[0,r] X [0,oo). Then we have 



4 



Lemma 2.2 The weak solution u{t, x) is obtained by getting it from u{t, x) of ([s]) by the trans- 
formation u{t, x) := u{t, X — f3{t)), where 

■ , , 1 du , 
/3(t) = -^(t.0+). 

This technique is also used in Section |4j 

2.3 Existence and Uniqueness of the Truncated Solution 

First, in order to control the nonlinear drift term we shall truncate the i^^-norm of the solution 
(which is shown as equivalent to truncating the nonlinear term, or P{t)), and work with the 
truncated solution u^{t,x), namely, the solution to 

u^{t,x)= p{t,x,y)uo{y)dy 
Jo 

+ / / $is)q{t-s,x,y)^L{\\u^is,-)\\H)u^{s,y)dyds (4) 
Jo Jo 

rt roo 



/ / pit - s,x,y)u^{s,y)d^tix)ds, 
Jo Jo 



+ 

lo Jo 

where '■ [0,oo) — t- [0, 1] is as a smooth monotone decreasing function that satisfies X[o,L] ^ 
< X[0,L+1]- 

The existence and uniqueness of the truncated solution is proved by a Picard-type iteration 
on i/^-spaced combined with similar calculations in Lemma 3.3 of [7]. Unlike in Section [s] and 
Section [4j in this problem we need not worry about the existence of the spatial derivatives, or 
in particular, 

du , 

which is the right hand side in the Stefan boundary condition, because the stochastic perturba- 
tion is smooth in space. By calculations of such an iteration and the structural results about 
i^'^-space (see Section 5.1 of [6]), combined with Lemma 3.3 of [7], we have 

Lemma 2.3 Fix L > 0. Then the solution u^{t,x) to Q exists and is unique. 

2.4 Relaxation of the Truncation 

Define the stopping time 



Also, define 



and 



:=inf{t > : ||u^(t,-)lk > L} 



T := lim T^, 

L— >oo 



u{t, x) := lim u {t, x). 

L—^oo 



We then have 
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Lemma 2.4 



and 



limllnft, AWh = oo 



lim 



du 
dx 



it,o+) 



oo. 



Finally, we have 

Lemma 2.5 The solution u{t, x) to ([s]) exists and is unique. 

Combining all the lemmas, we showed the main theorem. Note that the idea of first stopping 
\f3{t)\ from growing too large (that is, exceeding a fixed L), then using this to control the 
nonlinear drift term, and finally relaxing this truncation and obtaining a global stopping time 
T is important, and is also used in Section |4] when we study a stochastic Stefan problem driven 
by a scaled space-time Brownian noise. 



3 Boundary Regularity of the Stochastic Heat Equation 

In the previous section we have shown the existence and uniqueness of a stochastic Stefan 
problem with a spatially-colored and Brownian-in-time noise. Since we would further study a 
stochastic Stefan problem with space-time Brownian noise under certain regularity conditions, 
it is necessary that we first understand the effect of such noise on the regularity of the boundary, 
namely, the differentiability of the solution at the moving boundary, in addition to the existence 
and uniqueness of the solution itself. This task is critical because the Stefan boundary condition 
of such a problem involves the spatial derivative of the solution at the boundary, and since the 
noise is Brownian in space (as well as in time), the solution in general may not have a spatial 
derivative everywhere except at the boundary. 

Therefore, to simplify the problem, we first in this section consider a stochastic heat equation 
of u driven by a multiplicative space-time Brownian noise, so that the noise vanishes at the 
boundary (where u = 0), and we would expect that u is differentiable just at the boundary. 
Specifically, by removing the shift effect of the moving boundary and studying a stochastic heat 
equation of this kind, we look to understand 

(1) under what sense (or, in what normed space) the solution exists, and the connection 
between such a norm or space and the differentiability of the solution at the moving 
boundary; 

(2) in what sense (P-a.s.? in L^? etc.) the Stefan boundary condition holds; 

(3) the spatial regularity of the solution (Holder continuity? with what parameters?) which 
may guide us on the study of a Stefan problem with a space-time Brownian noise in the 
next section, in particular, the effect of a boundary shift on the iteration of the Ito integral. 
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3.1 Problem and Main Theorem 

Fix a probability space {0,, 3", P), and suppose VF:ilxMxR— )-Risa standard 2-dimensional 
Brownian sheet. Consider the (formal) stochastic heat equation of u(t,x) with a multiplicative 
space-time Brownian noise on [0, T] x [0,oo), under a Dirichlet boundary condition: 

du d'^u d^W 

m(0, x) = uo{x), 

u{t,0) = 0,Vx < 0,t G [o,r]. 

From the classic work of [7] by Walsh, the weak solution to the formal equation ([5| is 
equivalent to the evolution equation 



u{t,x)= p{t,x,y)uoiy)dy + / p{t - s , x , y)u{s , y)W {dyds) 
Jo Jo Jo 

where t G [0,T],x > and p{t, x, y) is the corresponding kernel as defined in the previous section. 
Then we have the following theorem as the main conclusion of this section: 

Theorem 3.1 (1) The solution u{t,x) to ([s]) exists and is unique with respect to a normed 
space; 

(2) P-a.s., for all t G [0, T], u{t,x) is differentiahle at x = 0, and 

du uit x'] dv /■* dv 

— (t,0) = lim = / ^(t,0,y)uo{y)dy+ / ^{t,0,y)u{s,y)W{dydsy, 

ox x\o x Jq dx Jq Jq dx 

(3) For t £ [0,T], define v{t,x) := u{t,x)/x for x > and v{t,0) := limx\ov{t,x), then 
F-a.s., v{t,x) is (I — e, g — e) -Holder continuous on [0,T] x [0, 1] for e > 0. 



Theorem |3.1| is proved in Section 3.3 using the integral regularities of a newly defined ker- 
nel p{t, X, y) in Section 3.2 combined with a newly defined norm and an argument based on 



Kolmogorov's Continuity Theorem (see [3]) and Lemma 3.3 of [7J in Section 3.3 



3.2 Integral Regularities of Kernel p(t,x,y) 



The existence, uniqueness, and regularity of the solution described in Theorem 3.1 are based on 
a number of integral regularities of a newly defined kernel p{t, x, y). We present and prove them 
in this section. 

Define a new kernel p{t, x, y) as 

y 

p{t,x,y) := -p{t,x,y),\/x > 

X 

and 

dp 

P{t,0,y) ■■= lim. p{t,x,y) = y — {t,0,y), 

x\0 OX 

then we have 
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Lemma 3.2 (l)'ix> 0, 



(2) Vx,ye [o,i],fe [o,r], 



/>oo 

/ f{s,x,y)dy< 
Jo 



C_ 



n\p{s, X, z) - p{s, y, z)f dzds < Ct{x -y)^; 
_ 



(3) Vx > 0,s,t G [0,T], 

/ / f{r,x,y)dydr+ I / \p{r + {t - s),x,y) - p{r,x,y)f dydr < DT\t - s\'i . 

Js Jo Jo Jo 

Proof Wc only need to prove the above facts for [0, T] x (0, 1], since by Fatou's Lemma they 
can be extended to the cases for x,y = 0. 

Throughout the calculations we repeatedly use the following facts: 



(a) 



(b) if / and g are even, then 

f'OO -j^ POO 

/ g{y)[f{y-x) + f{y + x)\dy= - I g{y)[f{y-x) + f{y + x)]dy 

Jo ^ J-oo 

= o/ [9{y-x) + g{y + x)]J{y)dy= I [g{y - x) + g{y + x)]f{y)dy; 

^ J-oo Jo 

(c) using the fact that 1 — exp(— x) < 1 A x, we have for ^4 > 0, 



/' 

Jo 



l-exp(-f) 



ds , Ads 
ds < —i=ds + / — 1= 
s Jo Vs Ja 



4VA. 



Then 

(1) 



f 

Jo 



p^{s,x,y)dy 



Ci 



2 / y 

x-^s Jo 



+ e 



2e 



dy 



Ci 



^/ 9 9n —y ^ 9 ^ +y 



2Ci _yl , 2Ci / _2£ 

e s dy-\ K- 1 — e » 

x'^s \ 

^ 2CirJ|) ^ 2Cir(|) _3Cir(^) 



dy 



f 

Jo 



CX) 2 

ye a dy 
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(2) Suppose X <y, and define h := y — x. Then from above, 



roo CiT ( 

/ \P{s,x,z)-p{s,y,z)fdz= y= 

Jo V'S 



4 + s 



1 — e s 1 — e" 
+ 



x"^ {x + hY 



Ix 



where 



2Ci 



x{x + h)s 

2Ci 
x{x + h)s Jq 

AC 



e 2s — e 2s ] [ e 2s — e 2s ] dz 



z 



^2 f _(x+_^-zf_ (x+^+zf 

s -)- e 3 



e 4s e 



1 J _/t- 

e 4s 



2 /"OO 





x{x + h)s 
x{x + h)'^/s 

Cir(i) _^ f 1 



2 + X + 



e " dz — e " / 
7o 



z^ + 



e o dz 



e 4s 



s + 2{x + 



s + 2 



e o 



I x{x + /l) 

Note that from (c) above, we get 



■s/s x{x + h) 



dz 



f- 

Jo 



e 4s 



ds < 2h. 



Then we need to consider 



J{t,x,h) J \ ^ 



\ — e a 1 — e " — e 

H — -^T^^ 2e~4^ 



x^ {x + K)"^ 



x{x + K) 



ds. 



Now we use two methods to bound J{t,x, h). 

2 

(I) Using the fact that x — \ <1 — exp(— x) < x, we get 



Jit, X, M < / ^ (l - e-^) ds + + h) '-^ds <4h + T (1) 
(II) Define a function 



^t,x) :-- 



t 1 — exp ( — 



-ds 



Then 



dx 



it,x] 



l-e~^ 2 
2t ^ < -. 
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Then we have 



J{t, X, h) < 



<^{t, x) + ^{t, x + h) -2^ [t,x + 



h 2 

< - + 



h 



h'^y/i 2h 

+ . , . < — + 



+ 



x{x + h) Jq 



s I 1 — e 4s I 



2 \x 2; + I y x{x + h) X x{x + h) 

Combining (I) and (II), we have that for some Ct > 0, when x < h?/'^, using (I) and 
we get [p(s, X, z) — y, z)]^ < Crh^^^', when x > /i^/^, using (II) and we get 
/o°° [p(s, X, z) - p{s, y, z)f dz < Cxh^/'^. 

(3) Suppose s < t, and define k := t — s. The first integral on the left hand side is bounded 
by C/Vk using (1). Now, from (1), 



j j \p{r + k, X, y) - p{r, x, y)f dy = SCiT s + k 

where (defining k' := 5^^^) 



k + Js 



2Ci 



x'^^r{r + k) Jo 



2Ci 



x'^y^r{r + k) Jo 
4Ci 



x'^y^r{r + k) 
1' 



g 2(r+fc) _ g 2{r + fe) 
g r + fc' -|- g r + fc' 



00 y2 

{x^ + y'^)e~^ dy - e ''+2 



g 2r 



e 



2r + fc 
r + fc' 



dydr 

(«+2?Tfc") 
-|- g r + k' 



dydr 



= 2CiT 







A:2 



y/r{r + A;) 

1 r(r + A;) 



(2r + fc)2 
1 



(2r + A:)2 



1 



+ ) 6 dy 



dr 



Therefore we have for some Dt > 

fS />oo 



^dr > SCiT 1 „ I / 3 ■ 

fcA 2 \^ J Jo {r + k)2 



Vr + k' + 
rdr 



-{r + A;')2 



/ / [v{f + k,x^y) — p{r^x,y)] dy < D^k'^ JOi 
Jo Jo 



3.3 Proof of the Main Theorem 



Theorem 3.1 is proved in two steps. First we show that the solution to ^ exists and is unique 
in a normed space, where the norm is defined so that in the second step the regularity results of 
the solution can be derived by using the defined norm and the integral regularities of the kernel 



p{s, x, y) in Lemma 3.2 In other words, the norm we defined in the next part characterizes the 
essential component from which we derive the desired Holder continuity of the solution and its 
differentiability at the boundary, which are shown by using Kolmogorov's Continuity Theorem 
in the second step. 
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3.3.1 Existence and Uniqueness of the Solution 



Existence and uniqueness of the solution is shown by a Picard-type iteration. Consider the 
iteration 



Un+l{t,x) 

Or, if we define 



p{t,x,y)uo{y)dy + 



t poo 



JO 



Vn{t,x) := ,vo{x) : 



p{t, X, y)un{s, y)W{dyds). 
uo{x) 



(6) 



X 



then 



Vn+iit,x)= p{t,x,y)vo{y)dy + / p{t,x,y)vn{s,y)W{dyds). 



t noo 



JO 



(7) 



Fix p > 1. Suppose {/(x)}^>o is a stochastic process. Define a norm 



2p:=supE[/2P( 

x>0 



X 



Then we have 



Lemma 3.3 For all p > 1, the solution v{t,x) to ([T]) exists and is unique in the space defined 
by the norm \\ ■ \\2p. 



Proof We have 

Un+lit,x) - Un{t,x) 
X 

or equivalently, 

Vn+lit,x) - Vn{t,x) 



t j-OO 

JO 



t fOO 



Un{s,y) - Un-lis,y) 



p{t — s, X, y)W{dyds) 



JO 



[vn{s,y) - Vn-iis,y)]p{t - s,x,y)W{dyds). 



Define Hn{t) := ||fn(i, •) — Vn-i{t, Olbp- Then we have from Lemma 3.2 (1) that 



i^n+l(t) <supCpE 

x>0 

< sup Cp 

x>0 



Jo 

t /"OO 



<c'pt- 



Jo 

t 



{vnis,y) - Vn-i{s,y)f p^{t - s,x,y)dyds 

Jo JO 



p^{t — s, X, y)dyds 



E 



{vn{s, y) - Vn-i{s, y)f'^ p^{t - s, X, y)dyds 



ds < CpT 



Hjs) 







ds. 



where the first inequahty comes from Burkholder inequahty, the second inequahty comes from 
Jensen's inequahty that for p > 1, 



p-i 



aPb, 
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and the third inequahty comes from Lemma |3.2| (1). By Lemma 3.3 of [7], 

^Hn{t) < oo 

n 

and the convergence is uniform on compacts. That is, if < t < T, then Y2n ^nit) < Cp^T < oo 
where Cp^T is a constant dependent only on T. This gives immediately by Picard-type iteration 
that the solution v{t,x) exists and is unique with respect to || • \\2p- Moreover, 

\\v{t,-)hp < Wvohp + '^Hnit) < Dp^T < oo.D 

n 

Define u{t,x) := xv{t,x), then u{t,x) is the unique solution to 

3.3.2 Regularity of the Solution and Differentiability at the Boundary 

In this part we shall prove the differentiability of the solution at the boundary by giving a P-a.s. 
limit or the spatial derivative at the boundary. This is shown by using Kolmogorov's Continuity 
Theorem combined with the integral regularities shown in the previous parts and the existence 
of the solution under the || • \\2p norm. In fact, a stronger result is shown, namely, the solution 
is P-a.s. Holder continuous with parameters (| — e, g — e) on [0, T] x [0, 1] for e > 0, which can 
also be used to evaluate the impact of its spatial regularity on the boundary shift effect required 
in the stochastic Stefan problems drive by space-time Brownian noise. 

Lemma 3.4 P-a.s., for all t € [0,T] the solution v{t,x) to ([7]) is continuous at x = and 



du , . , , , 
— t,0 = hmv t,x 

ox x\0 



dp 
dx 



{t,0,y)uo{y)dy + 



t fOO 



JQ 



dp 
dx 



{t,0,y)u{s,y)W{dyds). 



Moreover, fort G [0,T] define v{t,0) := limx\ov{t,x), thenF-a.s., v{t,x) is (| — ejg— e)- 
Holder continuous on [0,T] x [0, 1] for e > 0. 

Proof We only need to consider the Brownian term. Fix p>l. For t G [0,T],x > 0, define 



h{t,x) :-- 



1 



t fOO 



JO 



uis, y)p{t - s, X, y)W{dyds) 



t l-OO 



JQ 



v{s, y)p{t - s, X, y)W{dyds) 



and 



/i(t,0) := 



t i-oo 



JQ 



u{s,y)^{t - s,0,y)W{dyds) 



t /"OO 



Jo 



v{s, y)p{t - s, 0, y)W{dyds). 



Then for t £ [0, T], < x, y < 1, we have 
E[ihit,x)-hit,y)fP] <CpE 



< Cr, 



JO 

t roo 



Jo 



v'^{s, z) [p{t — s, X, z) — p{t — s, y, z)] dzds 

-iP-i 



[p{t — s,x, z) — p{t — s,y, z)] dzds 

t /-oo 



JO 



E [i^^^(s, z)] [p{t — s, X, z) — p{t — s, y, z)]^ dzds 



t fOO 



< CpDp^T 

<C;r(x-y)t 



[p{t — s, X, z) — p{t — s, y, z)]^ dzds 



JO 

p 
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where the first inequahty comes from Burkholder inequality, the second comes from Jensen's 
inequality, the third comes from the previous part about uniform boundedness of •)||2p on 
[0,T], and the last one comes from Lemma 3.2 (2). 
Similarly, We also have for < s < t < T, x G [0, 1], 



E [{h{t, x) - h{s, x)fP] < 22p-i[/n(s, t, x) + /i2(s, t, x)] 



where 



Iii{s,t,x) :=E 
Ii2{s,t,x) :=E 



s />oo 



JO 

t roo 



s Jo 



v{r, y) [p{t - r, x, y) - qs-r{x, y)] W{dydr) 

2p 



2p 



v{r, y)p{t — r, x, y)W [dydr] 



by Jensen's inequality (^)^^ < ^^^^±^. Now, 



/ii(s,t,x) < CpE 



< a 



2 

rs roo 

/ / v^ir,y)\p{t - r,x, 
Jo JO 

[p{t - r, X, y) - p{s - r, x, y)f dydr 

oo 

E \v'^^{r, y)] [p{t — r, x, y) — p{s — r, x, y)]^ dydr 



Jo 



,y) -p{s - r,x,y)] dydr 



JO 



using Burkholder 's inequality, Jensen's inequality, and Lemma 3.2 (3). Similarly, using Burkholder 's 
inequality, Jensen's inequality, and Lemma 3.2 (1), we have 

/l2(s,i,x) < C2,p,T{t - 

Summing things up, for all p > 1, we have for < s,t < T,0 < x,y < 1, 

E[{h{t,x)-h{s,y)fP] <Mp,T(x-y)§ +A^p,T(t-s)2. 

By Kolmogorov Continuity Theorem, the above calculations imply that P almost surely, Ii{t, x) 
is (7, ;g)-H61der continuous on [0, T] x [0, 1] for all 7 < 1/4, /3 < 1/6. This implies that Ii{t, 0) = 
limx\o Ii{t,x), and that 0) is 7-Holder continuous on [0, T] for all 7 < 1/4. □ 



4 A Stochastic Stefan Problem with Space-Time Brownian Noise 

In the previous section we considered the stochastic heat equation of u driven by a multiplicative 
space-time Brownian noise and showed the existence and uniqueness of the solution under 

normed space defined by || • \\2p for all p > 1 and the P-a.s. differentiability at the boundary and as 
a by-product, the Holder continuity of u. The proof provides us with important hints on how to 
prove the existence and uniqueness of a stochastic Stefan problem driven by space-time Brownian 
noise, since we shall also prove that the Stefan boundary condition holds simultaneously at the 
moving boundary. 
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Fix a probability space ($7,9", P), and suppose Vl^rrixMxM— J-Mis the standard 2- 
dimensional Brownian sheet. Then the stochastic Stefan problem of u{t,x) has the general 
(formal) form 

52 



+ a{t, X, u{t, x))-Q^yx > Pit), 



du 

m 

u{t,x) = 0,Vx < /3(t), 
u(0,x) = uo{x), 

A/ N u(t, x) 

p${t) = hm -^-^ 

x\l3{t) X 



-a.s. 



for < t < r where r is some well defined stopping time. 

The introduction of a moving boundary brings about a number of challenges that do not 
exist in the stochastic heat equation problem, and one of them is to control the boundary shift 
effect, namely, if we iterate on the boundary /3n{t), then dn{t) := (3n+i{t) — fin{t) is the shift of 
the boundary between iterations. We look to control 



Jo Jo 



p{t - s,x,y + dn{s)) 



'0 Jo 
whose variance is 
rt 



a(s, y + dn{s),u{s, y + dnjs))) 
y + dnis) 



p{t - s,x,y) 



11 

Jo Jo 



~i. , , / ^^o■(s,y + (in(s),M(s,y + (i„(s))) a{s,y,u{s,y)) 

p{t -s,x,y + dn{s)) , , , . p{t - s, X, y) 

y + dn[s) y 



W{dyds) 



dyds. 



For multiplicative noise where a{t,x,u) = u, the variance above in terms of dn is completely 
determined by the spatial regularity of u or u := u/x. From the previous section we know that 
spatially v is continuous only at the boundary and has (| — e) -Holder continuity for e > on 
[0, 1]. To evaluate the impact of the spatial regularity on the above variance in terms of dn, we 
let a{t,x,u) := x, and a calculation with two methods similar to Lemma 3.2 (2) shows that 

ft f°° 2/2 

/ / [p{t — s,x,y + d) — p{t — s,x,y)] dyds < Cj-d^ . 
Jo Jo 

Although the tightness of the inequality is not justified, it still gives strong evidence that it would 
be difficult to have an iteration on d„, even if a is as smooth as cr = x. Therefore, instead of 
working with a multiplicative noise, we in this section consider a stochastic Stefan problem drive 
by scaled space-time Brownian noise (j(x) where cr(x) satisfies certain regularity conditions: 



Definition A function a : [0, 00) 



is a regular scaling function if cr{x) is Lipschitz in x and 



cr(x) ~ Ax°' at with a > |. Equivalently, |o"(x)| < min{Ax", Bx} for some A, B > 0. 

Note that Ax"' provides enough smoothness to control the boundary shift effect as iterating 
on f3n{t) below and Bx provides enough control as x — t- 00 for other parts so that the noise does 
not grow too large at infinity. 
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4.1 Problem and Main Theorem 



Suppose a : [0, oo) 
problem of u{t,x) 



is a regular scaling function. Then we consider the stochastic Stefan 



92 



du 
dt 

u{t,x) = CVx < 
n(0,x) = uo(x), 



(8) 



p/3(t) 



lim 

x\l3(t) 



-a.s. 



for < t < r for some well defined stopping time r. 

As before we work on a weak formulation of ([s]) and its transformed equivalent evolution 
equation. Unlike in the colored noise case, we do not directly have the differentiability of the 
solution at the boundary, so in the iteration of the boundary j3 we cannot simply let it be the 
spatial derivative of the solution at the boundary, but the one similar to what was obtained in 
the stochastic heat equation, and finally we have an extra step to show that the two coincide 
P-a.s., or equivalently, the Stefan boundary condition in ([s]) holds. 

Theorem 4.1 The solution u{t,x), p{t) to ^ exists and 

where := inf{i G M+ U {0} : |/3(t)| > L}. Moreover, F-a.s., j3 G C'^{[0,t)) and ${t) is 
Q — -Holder continuous for e > 0. 



Theorem 4.1 is proved in Section 4.3 using the integral regularities of the new kernels p{t, x, y) 
and X, y) proved in Section 3.2 and Section 4.2 , The regularity result on the moving boundary 
is proved via using Kolmorogov's Continuity Theorem similarly to the proof given in Section[3.3[ 



4.2 Integral Regularities of Kernels p{t, x,y), q{t, x,y) 

As in the previous section, the proof to the main theorem is based on the integral regularities of 
newly defined kernels. Define two new kernels (where p is the same as in the previous section) 



y dp 
pit,x,y) := -p{t,x,y),p{t,0,y) := lim p{t,x,y) = y^{t,0,y) 

X x\0 ox 



q{t,x,y) := -^{t, x,y), q{t,0,y) 
X ox 



d^p 

lmiq{t,x,y) = y—^{t,0,y) 



y 

tVi 
= 



exp 



'2t 



Then we have the following integral regularity results about q{t,x,y), similar to those for 



p{t, X, y) as in Lemma 3.2 



Lemma 4.2 There exists C > 0, K > such that 
(1) 



\q{s,x,y)\dy < 



C_ 
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(2) forx,x + h£ [0,1], 



(3) for k>0, 



/ / \q{t - s,x + h,y) - q{t - s,x,y)\dyds < Khe; 
Jo Jo 

/ / \q{t + k — s,x,y) — q{t — s,x,y)\dyds < Kki. 
Jo Jo 



Proof First we observe that 

,2 



q{t,x,y) = —p+{t,x,y) - -p{t,x,y) = —p+{t,x,y) - -p{t,x,y) 
xs s xs s 



where 



P+it,x,y) :-- 



Ci 



e 2a -\- e 2s 



Therefore the calculations to prove the above facts are fairly similar to those in Lemma 3.2 
hence we only provide sketch here. 



(1) this can be shown by using facts (a) and (b) in Lemma 3.2, where 

g{y) ■■= \y\, 



f{y) ■■= |y|exp 

followed by a direct calculation of the integral after applying (b) with /, g. 



'2s 



(2) Similar to Lemma 3.2 (2), we still use a two-method approach to estimate the first part 
containing p+ : 



JO 



p+{t - s,x + h,y) _ p+{t - s,x,y) 
X + h x 



dyds < Chs . 



From Lemma |3.2| (1) and (c) we also have 

/■* ds , 
/ — / \(x + h)p(t — s,x + h,y) — xp(t — s,x,y)\dy < C'h. 
Jo s Jo 

Combining the above two facts, we have there exists K > such that 

/ / \q{t — s,x + h,y) — q{t — s,x,y)\dyds < Khe . 
Jo Jo 

(3) Similar to Lemma 3.2 (3), we first estimate the first part containing 

[p+{t + k — s,x,y) — p+{t — s, X, y)] dyds < Ck^ . 



\p{t + k — s, X + h,y) — p{t — s, X, y)\dy < C'k^ . 



From Lemma |3.2| (1) and (c) we also have 

'•^ xds 
s 

Combining the above two facts, we have there exists K > such that 











/ / \q{t + k — s,x,y) — q{t — s,x,y)\dyds < Kki.n 
Jo Jo 
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4.3 Proof of the Main Theorem 

As in the previous sections, we work with the weak formulation of the solution and its equivalent 
evolution equation. Define u{t,x) := u(t,x + /3{t)), Then the equivalent evolution equation gives 

/•oo 

u{t,x)= p{t,x,y)uo{y)dy 
Jo 

/"CxD 

+ 11 Ks)qit - s,x,y)u{s,y)dyds (9) 



Jo 

t /-oo 



JO 



p{t - s,x,y)a(y)W(s{dyds), 



where Wj^ is the Brownian sheet obtained by shifting W spatially by x — )■ x + f3{t). The main 
theorem is proved by adopting the following strategy: 

(1) prove the existence and uniqueness of the truncated solution u^{t, x) and its corresponding 
/3{t) in II • ||2p for < t < for a fixed L > by an argument based on Picard-type 
iterations; 

(2) based on (1), show that P-a.s., the Stefan boundary condition 

hm — = pP(t) 

x\0 X 

holds by using Kolmogorov's Continuity Theorem on the drift term and the Brownian 
term; 

(3) let L — 7- oo and obtain the solution u{t,x); transform u back to u and obtain a weak 
solution of ([s]). 

4.3.1 Existence and Uniqueness of the Truncated Solution 



In this section we first show the existence and uniqueness of (3 that satisfies (10), and then 
truncate f3{t) by a fixed L > so that the drift term of the solution is controlled from growing 
too large. Then we show the existence and uniqueness of the solution under the truncation. 
Defining Wjj as before, then we have 



Pf3{t) = 1^ '^{t,0,y)uo{y)dy 

+ r r ^{t-s,0,y)a{y)Wp{dyds), 



JO 



(10) 



Lemma 4.3 There exists a unique (3{t) that satisfies (10). 
Proof Consider the following iteration on /3(t): 

Pl^nit) = I -^{t,0,y)uo{y)dy 







+ [ r p{t-s,0,y)^Wp^_,{dyds) 

Jo Jo y 



17 



Fix p > I and define 



Also define 



Hn{t) :=E (/3„+i(t)-/3„(t))2p 



g{s,y) := p{s,0,y)y°' \ 
Then letting dn{t) := /3n(i) — f3n-i{t), we have 

Hnit)<KiElj^ [^{t-s,y + \dn{s)\)-g{t-s,y)fdyds\ 



< KiE 

< K2E 



\dn{s) 
1 rt 



JO 



dg 



(t- s,y + X\dn{s)\)dX 



\dn{s) 

1 rt 



dy 
(y + A|d„(s)|)2° 



dyds 



exp 



(y + A|d„(s)|)2 



|(i„(s) 



JO 



.5 exp 



dydsdX 



(t-s) 



it-s) 

- {y + X\dnis)\f 

(t-s) 



dydsdX 



\dnis)\' 



-ds )■ < Kr^E 



K.E 



(t-s)2-° J 
ir 



JO 



(t-s)^ 



drds 



dsdr 



Note that for the above argument to work we must have a > | , otherwise the first integral with 
respect to r is 00. Since | — a > — 1, we get that 



Also, define 



P/3o(i) 



dp- 



dx 

t roo 



{t,0,y)uQ{y)dy 



+ 



Jo 



then from Lemma |3.2| (1) we have 



p{t-s,0,y)^W{dyds), 

y 



E|/3o(t)pP < 00 



Therefore, by Walsh's Lemma 3.3 in [7J, we obtain that f3{t) as the limit of /3n(i) exists and 
is unique, and we also have 

E\${t)\^P < 00 

and its bound is uniform (that is, not dependent of t). □ 



Let /3(t) be the solution to (jlOj). Fix L > 0, and define '^l '■ [0,oo) — )• [0,1] as a smooth 
monotone decreasing function that satisfies X[o,L] < < X[o,L+i]- Now consider the following 
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iteration of u^{t, x): 



Un{t,^)= P-{t,x,y)uo{y)dy 
Jo 

rt roc Qp 



I3{s)'^{t - s,x,y)u^_^{s,y)^L{\$is)\)dyds 



(11) 



'0 ^0 

ft f'OO 

+ 11 P-{t - s,x,y)a{y)Wi3idyds), 
Jo Jo 



where the initial value 



fOO POO 

UQ{t,x):= p-{t,x,y)uo{y)dy + / p-{t - s, x,y)a{y)W^{dyds). 
Jo Jo Jo 



Then we have 



Lemma 4.4 The solution u^(t,x) to (11) exists and is unique. 



Proof By adopting a similar approach to the proof of Theorem 3.1, we first compute 



X 



t fCO 



JO 



P{s)q{t - s,x,y) 



Un{s,y) -u^_i{s,y) 



^Lms)\)dyds. 



Define 



Jn{t) '■= SUpE 

a;>0 



X 



2p 



From Lemma 4.2 (2), still use Burkholder's inequality and Jensen's inequality and we get 

s 2p 



ipt poo 
j J q(t-s,x,y) 

Jo Vi- s 



<{s,y) -M„_i(s,y) 



y 



dyds 



Also, from Lemma 3.2 (1) and < By, by following the same calculations as in the 

2p1 



previous sections we get 



supE 

x>0 



u^{t,x) 



< oo 



Therefore, by Walsh's Lemma 3.3 in [7J, we obtain that u^{t, x) as the limit of u^{t, x) exists 
and is unique, and we also have 



supE 

a;>0 



u^{t,x) 



X 



2p 



< OO 



and its bound is uniform (that is, not dependent of t). □ 
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4.3.2 The Stefan Boundary Condition Holds 



Lemma 4.5 Let and (3 be the unique solutions to (11) and (10). Then P-a.s., the Stefan 
boundary condition holds, namely, 



u^{t,x) 



lim 



a;\0 X 

and (5{t) is — e) -Holder continuous. 



P/3(i),V0 <t<T^ 



Proof This lemma is shown by using Kolmogorov Continuity Theorem and adopting the similar 
calculations as in the previous section. In fact, define for x > v{t,x) := u^{t,x)/x and 
v{t,0) := p${t), then we claim that P-a.s., v{t,x) is (| — e, g — e)-Holder continuous for fixed 
r > on [0, T] X [0, 1] with e > 0. 
Indeed, by Burkholder's inequality. 



E [\v{t, x + h)- v{t, x)\'^P] < KiE 



t l-OO 



JO 



-,2p 



[q{t - s,x + h,y) - q{t - s, x, y)]v{s, y)dyds 



+ K2E 



t roo 



JO 



[p{t - s,x + h,y) - p{t - s, X, y)] v (s, y)dyds 



By Lemma |3.2| (2) and Lemma 4.2 (2) about p and q, we have that by using Jensen's 



inequality as in the previous section, for < x < 1, 

E [\v(t,x + h) -v{t,x)\'^P] < Kha 



where K does not depend on t. Also, by Lemma 3.2 (3) and Lemma 4.2 (3), we have that still 
by using Jensen's inequality, for < t < T, 

E[\v{t + k,x) -v{t,x)\'^P] <Kk2. 

By Kolmogorov's Continuity Theorem, we have that for all 5 < 1/4, 7 < 1/6, v{t, x) is P-a.s. 
uniformly ((5, 7)-H61der continuous, which implies that P-a.s., 

limv{t,x) = p${t),yt G [o,r], 

x\0 

and I3{t) = v{t,0)/p is (| - e) -Holder continuous. 

Now, when < 00, we simply let T := r^. When = 00, we choose T := n for all 
n G N, and since the above statement holds for t G [0, n] for all n G N we have that it holds for 
t G [0,00). □ 



4.3.3 Relaxation of the Truncation 

Lemma 4.6 Define r := lim^^oo and u{t, x) := lim^^^oo u^ii, x). Then u{t, x) is the unique 
solution to (Isl). 
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Proof Fix t E [0, r) and define Lt := sup{|/3(r)| : < r < t}. Then we have Lt < oo. Therefore 
for < r < t we have that u{r,x) also satisfies the equation of n^*(r, x). By the uniqueness of 
P and M^*, this implies that for < r < i, 

u(r,x) = u^^(r,x) and lim ^ — = p$(r). 

x\0 X 

Therefore, for < t < r, tt(t, x) satisfies ([s]) and also 

U(t, X) A, X „ 

lim ^ ' ^ = pB(t)n 

x\0 X 

4.4 Numerical Simulation 

Since the existence and uniqueness of the problem is shown we can indeed simulate the solution 
and its boundary numerically. The simulation uses finite-difference Euler approximation scheme 
described in [1]. To guarantee numerical robustness in our simulation we assume the space and 
time increment steps satisfy At < (Ax)^/2. 

In this simulation we apply the following simulation parameters: 



p = -0.2 

uo{x) 



X + x"^ 



1 + a;4/16' 

aix) = — — . 
^ ' 1 + 4x 

The consequent numerical results are illustrated as follows. 

(1) Figure [2] illustrates the weak solution u{t,x)\ 

(2) Figure [3] illustrates the boundary derivative 

(3) Figure [4] illustrates the typical shape of the solution u{t, x) at a particular time t > 0, from 
which we see that u is smoother as x gets closer to the boundary (shifted and denoted by 
0), and is differentiable at the boundary. 



5 Summary 

In this section we summarize the main results obtained in the previous sections. Throughout 
this section we fix a probability space (r2,3", P), and suppose VFiJlxlRxIR— )-]Risa standard 
2-dimensional Brownian sheet. In this paper we studied 3 types of stochastic equations. 
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Figure 4: A typical shape of the solution u at some time t > 0. is smoother as x gets closer 
to the boundary (shifted and denoted by 0), and is differentiable at the boundary. 



5.1 A Stochastic Stefan Problem with Spatiahy Colored Noise 

The solution u{t, x) to the stochastic Stefan equation 

du d'^u 

u{t,x) = 0,Vx < p{t), 
n(0,x) = no(x), 

u{t, x) 



+ u{t,x)dS,tix),yx > P{t), 



fl2) 



pf3{t) 



lim 



-a.s. 



x\l3{t) X 

exists and is unique for < t < r := lim^^^oo where := inf{t G R^- U {0} : > i}. 

5.2 Boundary Regularity of the Stochastic Heat Equation 

The solution n(t, x) to the stochastic heat equation 

dx ^ ' dx 
■u(0,x) = uo(x), 

it(t,0) = 0,Vx < 0,t E [o,r] 



— + — ,Vx>0,tG [o,r], 



(13) 



exists and is unique; Moreover, define v{t,x) := u{t,x)/x, then P-a.s., 



^it,0,y)uoiy)dy + ^{t,Q,y)u{s,y)W{dyds); 
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If we further define v{t, 0) := lim^-x^o '^{t, x), then P-a.s., v{t, x) is (j — e, g — e)-Holder continu- 
ous on [0, T] X [0, 1] for e > 0. 

5.3 A Stochastic Stefan Problem with Space-Time Brownian Noise 

Let (T : [0, oo) — )• M be a regular scahng function (see the definition given in Section |4]). Then 
the solution u(t, x) to the stochastic Stefan equation 

du d'^u d'^W 

u(t,x) = 0,yx< Mt), 

~ (14) 
u(0,x) = uq{x), 

pm= lim P-a.s. 

x\l3(t) X 

exists and is unique for < t < r for < t < r := \\mL^ao where := inf{t G M-|- U {0} : 
\${t)\ > L}; Moreover, P-a.s., (3 G C^{[0,t)) and /3(t) is (| - e)-Holder continuous for e > 0. 
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